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Abstract 

The Large Magellanic Cloud, a satellite galaxy of the 
Milky Way, has heen observed with the High En¬ 
ergy Stereoscopic System (H.E.S.S.) above an energy 
of 100 billion electron volts for a deep exposure of 
210 hours. Three sources of different types were 
detected: the pulsar wind nebula of the most en¬ 
ergetic pulsar known N157B, the radio-loud super¬ 
nova remnant N132D and the largest non-thermal 
X-ray shell - the superbubble 30 Dor C. The unique 
object SN 1987A is, surprisingly, not detected, which 
constrains the theoretical framework of particle ac¬ 
celeration in very young supernova remnants. These 
detections reveal the most energetic tip of a 7 -ray 
source population in an external galaxy, and provide 
via 30 Dor C the unambiguous detection of 7 -ray 
emission from a superbubble. 

Introduction 

The origin of cosmic rays (CRs), the very high (VHE, 
> 10^^ eV), and ultra high (> 10^® eV) energy par¬ 
ticles that bombard Earth, has puzzled us for over a 
century. Much progress has been made during the last 
decade due to the advent of VHE 7 -ray telescopes. 
These telescopes detect ~ 10^^ — 10^^ eV 7 -rays from 
atomic nuclei (hadronic CRs) collisions with local gas, 
or from ultra-relativistic electrons/positrons (leptonic 
CRs), which produce 7 -ray emission by upscattering 
low-energy background photons (7). Indeed, a survey of 
the inner part of the Milky Way with H.E.S.S., an array 
of imaging atmospheric Cherenkov telescopes (2), re¬ 
vealed a population of supernova remnants (SNRs) and 
pulsar wind nebulae (PWNe) emitting 7 -rays with ener¬ 
gies in excess of 100 GeV (3). 

Here we report on VHE 7 -ray sources detected out¬ 
side the Milky Way, namely in the Large Magellanic 
Cloud (EMC). This satellite galaxy of the Milky Way 


has a stellar mass of about 4% of the Milky Way (6, 7). 
Located at a distance of « 50 kpc (8), it is an irregular 
galaxy seen almost face-on (9). Consequently, source 
confusion is much less of a problem than for the inner 
Milky Way, and there is less uncertainty in the distances 
of the sources. The EMC stands out among nearby 
galaxies for its high star formation rate per unit mass, 
which is about a factor of five higher than in the Milky 
Way (JO, 11), and contains the best example of a local 
starburst, the Tarantula Nebula. The EMC also harbors 
numerous massive stellar clusters and SNRs. Among 
the SNRs is a unique source, SN1987A, the remnant of 
the nearest supernova observed in modem times (72). 

High-energy 7 -ray emission from the EMC was de¬ 
tected by EGRET (13) and, more recently, by the Eermi 
Large Area Telescopes (EAT) (14), which revealed dif¬ 
fuse emission with an extension of several degrees in di¬ 
ameter, tracing massive starforming regions. VHE 7 -ray 
telescopes, like H.E.S.S., besides providing information 
on much higher energy CRs, have an angular resolution 
of a few arcminutes, which is substantially better than 
Eermi-LAT’s resolution at 7 -ray energies < 10 GeV. 
The good angular resolution allows H.E.S.S. to iden¬ 
tify individual sources in the EMC. As we will detail 
below, a deep H.E.S.S. observation revealed three lumi¬ 
nous sources in the EMC: the superbubble 30 Dor C, 
the energetic PWN N 157B, and the radio-loud SNR 
N 132D. Of these sources, only N 157B was detected 
previously in a 47-hours exposure (75). The observa¬ 
tions extend the scope of VHE 7 -ray astronomy by pro¬ 
viding examples of sources from a population outside 
the Milky Way. N 157B and N 132D belong to known 
7 -ray source classes, but both have distinguishing char¬ 
acteristics, N 157B being powered by the most energetic 
young pulsar, and N 132D being one of the oldest VHE 
7 -ray emitting SNRs. The superbubble 30 Dor C, how¬ 
ever, provides an unambiguous detection of a superbub¬ 
ble in VHE 7 -rays. Conspicuously absent from our list 
of three sources is SN1987A, despite predictions that it 


1 




Figure 1; Sky maps of the Large Magellanic Cloud, a) Optical image of the entire LMC (4). The boxes denote 
the regions of interest discussed in this paper. Colours denote levels of 3, 5, 10 and 20 a statistical significance of 
the 7 -ray signal, b) VHE 7 -ray emission in the region around N 157B. The green lines represent contours of 5, 10 
and 15 a statistical significance of the 7 -ray signal, c) XMM-Newton X-ray flux image of the region of 30 Dor C. 
The superimposed cyan lines represent contours of 68 %, 95% and 99% confidence level of the position of the 
7 -ray source. Diamonds denote the positions of the star clusters of the LH 90 association. See supplementary 
material for details on the X-ray analysis, d) VHE 7 -ray emission in the region around N 132D. The green lines 
represent contours of 3, 4 and 5 a statistical significance. The background of the 7 -ray emission (in panels b and 
d) was obtained using the ring background method (5). The resulting excess sky map is smoothed to the angular 
resolution of the instrument. 
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should be a bright 7 -ray source (76, 17). 

H.E.S.S. Observations 

We report on a deep, 210-hours H.E.S.S. exposure, tar¬ 
geted at the region of the Tarantula nebula — corre¬ 
sponding to 30 Doradus (30 Dor) — the largest star¬ 
forming region in the Local Group of galaxies. We re¬ 
constructed 7 -ray showers with an image-fitting analy¬ 
sis {18) and cross-checked with a multivariate analysis 
based on image parameterization (79, 20), with consis¬ 
tent results. In both analyses, a cut on the uncertainty of 
the reconstructed 7 -ray direction indicated an angular 
resolution of « 0.05°. 

Fig. [T^ shows an optical image of the EMC over¬ 
laid with TeV gamma-ray point-source significance con¬ 
tours. In this dataset, 613 7 rays with a statistical signif¬ 
icance of 33 a are detected from the PWN N 157B. Fig¬ 
ure provides a close-up view of the 7 -ray emission 
from N 157B. The diameter of N 157B of 100" (27) is 
of the order of the H.E.S.S. angular resolution. Further 
significant 7 -ray emission is detected to the South-West 
ofN157B. 

A likelihood fit of a model of two 7 -ray sources to 
the on-source and background sky maps establishes the 
detection of a second source at an angular distance of 9' 
(corresponding to 130 pc at a distance of 50kpc) from 
N 157B. The model consisting of two sources is pre¬ 
ferred by 8.8 tr over the model of one single source. 
Fig. & shows an X-ray image with overlaid contours 
of confidence of the source position. The position 
of the second source (RA = 5*'35™(55 ± 5)®, Dec = 
—69°11'(10 ± 20)", equinox J2000, 1 cr errors) coin¬ 
cides with the superbubble 30 Dor C, the first such 
source detected in VHE 7 -rays, and thus representing an 
additional source class in this energy regime. A 7 -ray 
signal around the energetic pulsar PSR J0540—6919 is 
not detected, despite the presence of an X-ray luminous 
PWN (22). A flux upper limit (99% confidence level) is 
derived at7^.^(> ITeV) < 4.8 x 10“^"^ phcm“^ s“^. 

Along with the clear detection of N157B and 
30 Dor C, evidence for VHE 7 -ray emission is observed 
from the prominent SNR N 132D (Fig. [^). The emis¬ 
sion peaks at a significance of about 5 a above a back¬ 
ground which is estimated from a ring around each sky 
bin. At the nominal position of the SNR 43 7 rays with 
a statistical significance of 4.7 a are recorded. 

The 7 -ray spectra of all three objects are well 
described by a power law in energy, ^{E) = 


(PN/{dEdtdA) = fl>o(F;/lTeV)”^ (Fig. ||i. The 
best-fit spectral indices and integral 7 -ray luminosities 
are summarized in Table [T] 

Even with a deep exposure of 210 hours, significant 
emission from SN 1987A is not detected, and we de¬ 
rive an upper limit on the integral 7 -ray flux of F.y(> 
ITeV) < 5.6 X 10“^^ phcm“^ at a 99% confi¬ 
dence level. 

Discussion of individual sources 

The three VHE emitters belong to different source 
classes and their energy output exceeds or at least equals 
that of their most powerful relatives in the Milky Way. 

30 Dor C 

The superbubble 30 Dor C stands out in X-rays as it con¬ 
tains, in the western part, an X-ray synchrotron-emitting 
shell with a radius of 47 pc, which makes it the largest 
known X-ray synchrotron shell (24-26). X-ray syn¬ 
chrotron emission, which indicates the presence of VHE 
electrons, is usually associated with 100 — 2000 year- 
old SNRs with radii smaller than 25 pc. In addition, 
the X-ray synchrotron luminosity of 30 Dor C is ten 
times that of the archetypal young SNR SN 1006 (24). 
The 30 Dor C shell also emits radio and optical radia¬ 
tion (27), and appears to have been produced by the stel¬ 
lar winds and supernovae in the OB association EH 90 
(NGC 2044) (28). 

The measured H.E.S.S. flux of 30 Dor C corre¬ 
sponds to a 1 —10 TeV 7 -ray luminosity of (0.9±0.2) x 
10 ^® erg s“^, with the best-fit position of the 7 -ray emis¬ 
sion lying in between the six identified sub-clusters (28). 
The TeV emission can be explained by the production 
of neutral pions due to collisions of hadronic CRs with 
the background plasma. Alternatively, the so-called lep- 
tonic emission scenario may apply, in which case the 
TeV emission is the result of Compton upscattering of 
low-energy photons to 7 -ray energies, by the same pop¬ 
ulation of electrons that is responsible for the X-ray syn¬ 
chrotron radiation (7). 

For the hadronic scenario, a combination of en¬ 
ergy in CRs (assumed to be protons) and den¬ 
sity of hydrogen atoms, nn, of Wpp = (0.7 — 

25) X 10®^ (nn/lcm”^)”^ erg is required (see SI.3). 
30 Dor C probably experienced ~5 supernova explo¬ 
sions (25), which likely provided ~5xl0^°erg in CR 
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Figure 2: Gamma ray spectra of N 157B, 30 Dor C and N 132D. Note that the spectral points of 30 Dor C are 
not corrected for spill-over emission from N 157B (see online supplement). Data points have 1 a error bars, upper 
limits are at the 99% conhdence level. The bottom panels show the residuals of the data points compared to the 
best-ht model. 


Source 

N157B 

30 Dor C 

N132D 

H.E.S.S. Identifier 

HESSJ0537-691 

HESSJ0535-691 

HESSJ0525-696 

Exposure Time 

181 h 

183 h 

148 h 

7 rays 

613 

74 

43 

Significance 

33.0 cr 

8 . 8 cr 

4.7cr 

Photon Index F 

2.8 ± 0.1 

2.6 ± 0.2 

2.4 ±0.3 

$(1 TeV) [ 10-12 cm-2 s"! TeV"!] 

1.3 ±0.1 

0.16 ±0.04 

0.13 ±0.05 

L^(l - 10 TeV) [10^5 ergs-i] 

6.8 ±0.3 

0.9 ±0.2 

0.9 ±0.2 


Table 1: Statistics and spectral parameters of the three sources. The exposure time is corrected for the acceptance 
differences due to different offsets from the camera centre. The signihcances of N 157B and N 132D are statistical 
signihcances of the 7 -ray emission obtained by using formula 17 of (23). The background was estimated from 
regions with similar offsets from the camera center as the on-source region. The signihcance of 30 Dor C is 
the significance by which a two-source morphology (N 157B and 30 Dor C) is preferred over a single-source 
morphology (N 157B only). The 7 -ray count and the flux of 30 Dor C are corrected for spill-over emission from 
N 157B (see online supplement). T is the photon index and $(1 TeV) the differential flux at 1 TeV of a power law 
fit to the energy spectrum. The luminosity is calculated for an assumed distance of 50kpc (S). The listed errors 
are statistical, 1 a errors. Systematic errors are estimated to be ±0.3 for F and ±30% for $(1 TeV) (75). 
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energy. Hence, the average gas density should be nn ^ 
20 ctn“^, which is higher than the density estimate of 
tt-h ~ 0.1 — 0.4 cm“^ based on the X-ray thermal emis¬ 
sion in the southwest (24, 29). However, locations of 
high densities may be present, if the X-ray thermal emis¬ 
sion comes from smaller radii than the dense outer shell, 
or if cool, dense, clumped gas survived inside the other¬ 
wise rarihed interior of the bubble (30). 

This hadronic scenario puts constraints on the CR 
diffusion coefficient, because the diffusion length scale 
should be smaller than the radius of the shell; Zdiff = 
\/2Dt < 47 pc for CRs around lOTeV. Therefore, 
D(10 TeV) < 3.3 x yr)-i cm^ s-\ which, 

given an age for the superbubble of a few million years, 
gives a much smaller diffusion coefficient than the typ¬ 
ical Galactic diffusion coefficient of 71(10 TeV) > 5 x 
10^®cm^s“^ (31). This small diffusion coefficient re¬ 
quires magnetic-field amplification combined with tur¬ 
bulent magnetic helds, as hypothesized by (32). 

X-ray synchrotron emission from 30 Dor C requires 
large shock velocities, Ushock ^ 3000kms“^ (55). As¬ 
suming that this shock originates from an explosion 
centered at the superbubble, we obtain a rough esti¬ 
mate of the age of the western X-ray shell of t = 
0.4i?/ushock ~ 6000 yr, assuming a Sedov expansion 
model (R = 2.8 x 10^(Et^/n-n)^^^ cm). Since the OB 
association is much older, this age most likely refers to 
a recent supernova explosion, whose remnant evolves 
in the rarihed medium of the superbubble. The Sedov 
expansion model then gives us a very low estimate for 
the density of nn ~ 5 x 10“^ cm“^ for an explosion 
energy of 75 = 10®^ erg. Although this is very low, it 
can occur under certain conditions (36). This model for 
the X-ray synchrotron shell can even be reconciled with 
the hadronic model of the TeV emission, if the rarihed 
medium also contains dense clumps. For the leptonic 
scenario for the TeV emission, the broad spectral energy 
distribution (SED, Fig. requires an energy in acceler¬ 
ated electrons of ~ 4 x 10"^® erg, and average magnetic 
held strength of 15 /iG, low compared to most young 
SNRs (57), but a factor three to four higher than the av¬ 
erage magnetic held in the EMC (38). 

Although at this stage we cannot rule out either the 
leptonic or the hadronic scenario, the H.E.S.S. observa¬ 
tions reveal that the conditions inside the superbubble 
must be extreme: the hadronic scenario requires loca¬ 
tions with high densities and a high degree of magnetic 
turbulence, whereas the leptonic scenario requires the 
stellar cluster to be extremely rarihed. Moreover, the 


7 -ray and X-ray observations suggest active particle ac¬ 
celeration by a very large, fast expanding shell. This 
may provide the right conditions for accelerating some 
protons to energies exceeding 3 x 10^® eV, which is the 
maximum energy detected for Galactic CRs. These ob¬ 
servations, therefore, lend support to the view expressed 
in (52, 59, 40) that superbubbles may provide the right 
conditions for particle acceleration to very high ener¬ 
gies, because they are thought to contain very turbu¬ 
lent magnetic helds and they are large enough to contain 
VHE particles for up to millions of year. 

In the Milky Way, the most closely related object to 
30 Dor C is the stellar cluster Westerlund 1 (47), which, 
however, has a completely different X-ray morphology. 
More importantly, it is not clear whether the 7 -rays orig¬ 
inate from the cluster wind itself, a PWN or from the nu¬ 
merous supernovae that exploded inside Westerlund 1 in 
the recent past. Since a large fraction of supernovae are 
thought to go off in superbubbles, this hrst unambigu¬ 
ous detection of VHE 7 -rays from a superbubble may 
have broad implications for the circumstances in which 
a large fraction of CRs are accelerated. 

N157B 

The source HESS J0537—691 is coincident with 
the PWN N157B, which sutTounds the pulsar 
PSRJ0537-6910. PWNe are nebulae of ultra- 
relativistic particles driven by highly-magnetized, 
fast-rotating neutron stars that convert a considerable 
amount of their spin-down energy into a particle 
wind. The archetypal Crab nebula is one of the 
brightest sources of non-thermal radiation in the sky 
and powered by the pulsar with the highest spin- 
down energy known in the Milky Way (42). With 
comparably extreme rotational energy loss rates, 
N 157B (75 = 4.9 x 10®®ergs“^) and the Crab nebula 
(75 = 4.6 X 10®®ergs“^) appear to be twins. The 
study of N 157B thus provides the unique opportunity 
to compare two extreme PWNe, and to disentangle 
object-specihc and generic properties. 

Given a population of ultra-relativistic electrons and 
positrons forming the PWN, the X-ray luminosity is 
determined by the strength of the magnetic field and 
the 7 -ray luminosity by the intensity of radiation helds 
which serve as targets for the inverse Compton upscat- 
tering. If the radiation helds are known, the magnetic 
held can be inferred from the combination of X-ray and 
7 -ray measurements. N 157B is likely associated with 
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Figure 3: Spectral energy distribution of 30 Dor C and N 132D. The 30 Dor C X-ray data are from (26). For 
N 132D, radio data is from (S3), and X-ray limits are from (34) and from re-analysed Chandra data. Both leptonic 
(dashed lines) and hadronic (solid lines) models are shown. For further details on Fermi-LAT data and spectral 
modelling, see S 1.2 and SI.3. 


the LH 99 star cluster (27, 43, 44), and therefore em¬ 
bedded in strong infrared radiation helds (see SI.5). In 
this environment, the magnetic held in the PWN must 
be rather weak, not larger than 45 pG, in order to ex¬ 
plain the multiwavelength data (Fig. |^. When con¬ 
sidering the region from which the hard X-ray emis¬ 
sion is coming, the total energy in the magnetic held 
is Wb, tot = 1-4 X 10'^^ erg - an order of magnitude 
smaller than the energy in >400 GeV electrons. The de¬ 
rived maximum magnetic held is also much lower than 
that inferred for the Crab nebula (^124/rG (45)), and 
suggests at least a factor ^7 lower magnetic pressure. 
As most of the electrons that radiate in the Chandra, X- 
ray and H.E.S.S., TeV domains have very short lifetimes 
(<300 years), the energy in ultra-relativistic particles in 
N 157B can be inferred independently of the spin-down 
evolution of the pulsar. For the model shown in Fig.|^ a 
constant fraction of 11 % of the current spin-down power 
of N 157B needs to be injected into the nebula in the 
form of relativistic electrons (compared with 50% for 
the Crab nebula under the same model). This fraction 
converted into X-ray and TeV emission is rather insen¬ 
sitive to the spectral index of injected electrons and the 
spin-evolution or braking index of the pulsar and only 
relies on the association of N 157B with LH 99 (see SI.5 
for more information). 

In this high-radiation field scenario, the situation 
for the Crab nebula is very different from N 157B. Not 
only is the best-ht electron spectrum of N 157B harder 
(Fe = 2.0 vs. 2.35), exhibiting a lower cut-off energy 
(Ec = 100 TeV vs. 3.5 PeV), but much of the spin- 


down energy of N 157B is also hidden and is not carried 
by ultra-relativistic particles or magnetic fields. The re¬ 
mainder of the available rotational energy is likely to be 
fed into electrons with energies <400 GeV that radiate 
at lower photon energies, adiabatic expansion, and/or 
particles escaping into the interstellar medium via diffu¬ 
sive escape (e.g., (46)). It therefore appears that N 157B 
is such a bright 7 -ray emitter because of the enhanced 
radiation fields, despite the fact that it is apparently a 
much less efficient particle accelerator than the Crab 
nebula. 

N132D 

In addition to the two unambiguously detected sources, 
we find strong evidence for a third source at the posi¬ 
tion of the core-collapse SNR N 132D. N 132D is a SNR 
with strong thermal X-ray emission, which has been 
used to estimate a pre-shock density of nn « 2.6 cm“^ 
(34), a high explosion energy of ~ 6 x 10®^ erg (34), and 
an age of ~ 6000 yr, based on a Sedov model. Such X- 
ray bright SNRs are predicted to be 7 -ray emitters (48). 
N 132D is also luminous in the radio (33) and infrared 
bands (49). N 132D is often compared to the brightest 
radio source Cas A, which, like N 132D, is an oxygen- 
rich SNR. N 132D has a higher infrared luminosity (49), 
but its radio luminosity is 50% that of Cas A. This is 
still remarkable given that N 132D has a kinematic age 
of 00 2500 yr (50), whereas Cas A is ~ 330 yr old and 
declines in luminosity by about 0.8% per year. The ra- 
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Figure 4: Intrinsic SED of N 157B (black) and the Crab nebula (grey). The model shown assumes the same 
injection parameters as derived for the Crab nebula (i?min = 400GeV, = 3.5 PeV, Fe = 2.35 (47)). The 
magnetic field required to explain the Chandra data of N 157B in the highest possible radiation fields is 45p,G. 
A significantly better ht to the Chandra data is obtained with Fg = 2.0, but requires a much lower cut-off of 
Ec < 100 TeV. 


dio properties have been used to infer a magnetic held 
strength of ^ 40 fj,G (33). The discrepancy between 
the age estimate based on the X-ray emission and the 
kinematic age may indicate that the supernova exploded 
within a bubble created by the progenitor star’s wind 
before encountering the high density material it now in¬ 
teracts with. 

The 7 -ray Hux measured by H.E.S.S. trans¬ 
lates to a 1 — 10 TeV 7 -ray luminosity of 
(0.9 ± 0.2) X 10^®(fi/50kpc)^ erg/s. Assuming 
that the 7 -ray emission is caused by neutral-pion 
production, this luminosity implies an energy of 
10®^(nH/lcm“^)“^((i/50kpc)“^ erg in relativistic 
protons. A hadronic origin of the 7 -ray emission, there¬ 
fore, implies either a large CR-energy fraction of 17% 
of the explosion energy, for an estimated post-shock 
density of nn ~ 10 cm“^ (34), or the gas density is 
higher than the x-ray-based estimates. The latter is 
plausible given that N 132D appears to interact with 
dense, shocked interstellar clouds, seen in the optical 
and the infrared bands (49). It is interesting to compare 
N 132D to the most luminous Galactic SNR detected 
at TeV energies, HESS J1640—465: both SNRs are 


believed to interact with a wind-blown cavity wall, 
to possibly have similar ages and sizes (50, 51), and 
to have transferred a large fraction of their explosion 
energies into CRs. 

The bright radio synchrotron luminosity of N 132D 
and the tentative claim of X-ray synchrotron emission 
from this source (52) also raises the possibility that the 
7 -ray emission is caused by inverse Compton scatter¬ 
ing of low-energy photons. In and around N 132D the 
radiation energy density is dominated by the bright in¬ 
frared flux from dust inside the SNR, and can be roughly 
estimated to be at least rtrad ~ 1-0 eV cm“^. This lep- 
tonic scenario requires that the average magnetic-field 
strength needs to be ~ 20 /iG, somewhat lower, but 
still consistent with the equipartition value (see SI.3). 
However, this leptonic scenario critically depends on 
whether the 4-6 keV X-ray continuum emission indeed 
contains a significant synchrotron component. 

Whatever the emission mechanism for the 7 -ray 
emission from N 132D, it is an exciting new 7 -ray- 
emitting SNR, because its age lies in the gap between 
young (< 2000 yr) TeV-emitting SNRs, and old (> 
10000 yr) TeV-quiet SNRs. The latter can be bright 
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pion-decay sources, but their spectra appear to be cut 
off above ~ 10 GeV. N 132D provides, therefore, an in¬ 
dication of how long SNRs contain CRs with energies 
in excess of 10^^ eV. 

SN1987A 

SN 1987A, the only naked-eye SN event since the Ke¬ 
pler SN (AD 1604), has been extensively observed at all 
wavelengths from the radio to the soft 7 -ray band, pro¬ 
viding invaluable insights into the evolution of a core¬ 
collapse SNR in its early stage (55). 

It has been suggested that even in the early stages of 
the SNR development, the shock wave, which is heat¬ 
ing the dense circumstellar medium (CSM) structured 
by stellar winds of the progenitor star, should have led 
to efficient acceleration of VHE nuclear CRs, accompa¬ 
nied by strong magnetic field amplification through CR- 
induced instabilities (16, 54). In collisions of the CRs 
with CSM particles, 7 -rays are produced. Estimates for 
the 7 -ray flux (16,17) strongly depend on the magnetic 
field topology and on the properties of the non-uniform 
CSM (55), making flux estimates uncertain by at least a 
factor of 2 (16). 

Based on a nonlinear kinetic theory of CR accel¬ 
eration, successfully applied to several young Galac¬ 
tic SNRs, the volume-integrated 7 -ray flux at TeV en¬ 
ergies, F.y{> ITeV), from SN 1987A was predicted 
to be rising in time, and to have reached a level of 
« 2.5 X 10“^^phcm“^ s“^ in the year 2010 (16). An 
analysis with different assumptions on CSM properties 
and a more phenomenological approach to CR accel¬ 
eration resulted in a predicted flux of ~ 8 x 10 “^^ 
phcm“^ s“^ in the year 2013 (17). The H.E.S.S. upper 
limit F.^(> ITeV) < 5 x lO-^^phcm-^ s"! at a 99% 
confidence level, obtained from observations made be¬ 
tween 2003 and 2012, being below the aforementioned 
predictions and a factor of 3 below similar estimates for 
the year 2005, therefore places constraints on the mod¬ 
els despite their uncertainties. 

The H.E.S.S. upper limit on the 7 -ray flux trans¬ 
lates into an upper limit for the 7 -ray luminosity of 
L.y(^> ITeV) < 2.2 x 10^^erg/s, which can be used 
to derive an approximate upper limit on the energy of 
the accelerated particles, Wpp, for a given average target 
density. Multi-wavelength studies of SN 1987A suggest 
that the shock at the current epoch has reached and is 
interacting with the so-called equatorial ring, for which 
gas densities ranging from 10^ cm“^ to 3 x 10"^ cm“^ 


have been found (56). Thus one finds a conservative up¬ 
per limit, Wpp < 1.4x 10^®/“^ erg, where 0 < / < 1 is 
the fraction of accelerated particles that are interacting 
with the dense regions. This upper limit on the energy 
of accelerated CR particles corresponds to 0.15/“^ % 
of the explosion energy of 10 ®^ erg. 

Assuming a spherically-symmetric distribution of 
accelerated particles, one can estimate / ^ 0.2 with 
the geometry of the equatorial ring found in (57). This 
translates to Wpp < 9 X 10^® erg, implying that less 
than 1 % of the explosion energy is carried by acceler¬ 
ated CR nuclei. This fraction is rather small compared 
to typical values of ~ 10% for young SNRs (of ages 
~ 1000 — 2000 years), but is not unreasonable for a 
very young object like SN 1987A. 

Summary 

With the deep H.E.S.S. observations of the EMC, we 
have detected three luminous examples of CR sources 
in an external galaxy. These sources detected in 7 rays 
include a superbubble and counterparts to the most lu¬ 
minous sources in the Milky Way. N 157B provides a 
counterpart to the Crab Nebula, but its electron acceler¬ 
ation efficiency is five times less than for the Crab neb¬ 
ula, and its magnetic field pressure is seven times less. 
N 132D has been long regarded an older version of the 
brightest Galactic radio SNR Cas A, and is one of the 
most radio-luminous SNRs known. N 132D is also re¬ 
markable in that it is one of the oldest VHE 7 -ray emit¬ 
ting SNRs. With the three detected sources, we increase 
our understanding of the variety of VHE 7 -ray sources, 
which will likely require observations with the future 
Cherenkov Telescope Array (58), which should be an 
order of magnitude more sensitive than H.E.S.S. 
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Supplementary Materials 

A Materials and Methods 

A.l Correction factor for 30 Dor C flux 

The superbubble 30 Dor C and the PWN N 157B have 
an angular distance of « 0.15°, while the angular res¬ 
olution ( 68 % containment radius of the point-spread- 
function (PSF)) of H.E.S.S. is « 0.05°. The flux of 7 - 
ray sources is calculated from the excess of 7 -like events 
in a pre-deflned on-region. As the reconstructed flux 
of 30 Dor C has been estimated within a test region of 
0.07° around the center of the X-ray superbubble, a part 
of the 7 -rays from N 157B lies in this on-region due to 
the tails of the PSF. These additional 7 -rays will result 
in a larger reconstructed flux for 30 Dor C, and need to 
be corrected for in order to estimate the “true” flux of 
30 Dor C. 

The H.E.S.S. PSF is approximately radially sym¬ 
metric. As the reconstructed extension of N 157B is 
compatible with that expected for a point source, the dis¬ 
tribution of 7 -rays from N 157B is expected to be radi¬ 
ally symmetric. Five circular spill-over test regions of a 
radius of 0.07° are placed at the same angular distance 
as between N 157B and 30 Dor C, but rotated around 
N 157B, as shown in Fig.[^ to estimate the 7 -ray excess 
contributed by N 157B in the 30 Dor C on-region. 

In total, 150 ± 25 7 -like excess (using the reflected 
background method (5)) have been recorded in the test 
regions, with an area 5 times larger than the on-region. 
The average 7 -ray excess contributed by N 157B in 
the 30 Dor C on-region is thus estimated to be (150 ± 
25)/5 = 30 ± 5. Excluding the region of highest or 
lowest excess does not significantly change this num¬ 
ber, demonstrating the robustness of this estimate with 
respect to the location of the spill-over test regions. The 
excess from 30 Dor C before correction is 104 ± 14, 
implying a flux contamination of 29% ± 6 %. The mea¬ 
sured flux was correspondingly scaled down by a factor 
0.71 ± 0.06. As the reconstructed 7 -ray spectral index 
of the superbubble and the PWN are consistent within 
errors and the energy dependence of the PSF is small 
in the energy range discussed here, the correction factor 
is considered to be independent of reconstructed 7 -ray 
energy. 

Another method to calculate the true excess from 
30 Dor C is to use the spill-over test regions as back¬ 
ground regions. In total, 183 and 544 7 -like events (in- 
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Figure 5: On (solid circle) and test (dashed circles) re¬ 
gions to estimate the flux from 30 Dor C on top of the 
smoothed 7 -ray excess map. The dash-dotted circle de¬ 
notes the region of interest of N 157B. The inset shows 
the simulation of a point-like source. 


eluding true 7 -rays and 7 -like CRs) have been recorded 
in the 30 Dor C on-region and the test regions, respec¬ 
tively. The excess from 30 Dor C calculated in this way 
is 74 ± 14, consistent with the value derived from the 
previous method. 

A.2 Fermi-LAT analysis 

A search for GeV 7 -rays from the sources detected 
by H.E.S.S. has been carried out using data from the 
space-born Fermi-FAT instrument. The FAT is a pair- 
conversion telescope, operating in the energy range be¬ 
tween ~ 20 MeV and ~ 300 GeV. A description of the 
instrument- and mission-related details is given in (59). 
Based on the data from an 11 month observation period, 
the detection of diffuse emission from the Farge Magel¬ 
lanic Cloud, best described by either a two-component 
Gaussian template or an HII template, has been reported 
(74). Emission from individual point-like sources above 
this diffuse emission has not been detected in this data¬ 
set. 

To derive the 7 -ray flux upper limits shown in Fig¬ 
ure 4 in the main article, more than 5 years of spacecraft 
data (2008-08-04 - 2013-10-21) have been analyzed us¬ 
ing the Fermi Science Tools (FST) package with ver- 
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sion v9r31plQ Only events between (10 — 300) GeV 
have been considered, which results in an excellent an¬ 
gular resolution of < 0 . 2 ° and reduces the contribution 
from the diffuse emission component. Events within a 
20 ° by 20 ° region and with zenith angles below 100 °, 
centered on RA 05‘’17“36.0" and Dec -69°01'48.0" 
(J2000) were analyzed with a binned likelihood analysis 
using the P7Clean_V6 instrument response functions. 

Sources have been modeled according to the Fermi 
2 -year catalog (60), and additional point-like sources 
at the positions of 30 Dor C, N 132D, N 157B, and 
SN 1987A have been added. Note that significant dif¬ 
fuse emission has been detected towards N 157B (74). 
As the origin of this emission cannot unambiguously 
be attributed to one of the H.E.S.S. sources, the narrow 
Gaussian component (G2 in (74)) was removed from the 
model. The normalization of all sources within 15° as 
well as that of all diffuse sources has been left free in 
the fit. 

The 7 -ray flux upper limits have been derived in 
three energy bands: (10 — 30) GeV, (30 — 100) GeV, 
and (100 — 300) GeV, assuming a power-law in energy 
with spectral indices as reconstructed for the H.E.S.S. 
sources (see Table 1). In the case of SN 1987A only flux 
upper limits could be derived in the TeV regime, there¬ 
fore a spectral index of E = 1.8, as motivated by (76), 
is assumed. 

A.3 SED modelling of 30 Dor C and 
N132D 

The leptonic and hadronic model curves shown in Fig¬ 
ure 3 in the main paper have been obtained using a 
model for the time-dependent injection and interaction 
of electrons and protons (e.g. (57)). 

For 30 Dor C electrons and protons are injected at 
a constant rate and cool via synchrotron and Inverse 
Compton processes, and proton-proton collisions, re¬ 
spectively. Given the rather low densities under con¬ 
sideration, we ignore Bremsstrahlung and Coulomb 
losses. The energy-dependent proton-proton cross- 
section is taken from (67). In a purely hadronic sce¬ 
nario, a total energy in interacting protons of Wpp = 
0.7 X 10®^ (jiYi/l cm“^)“^ erg is required to explain the 
H.E.S.S. data (correcting for a 30% contamination from 
N 157B, and assuming a cutoff energy of 100 TeV and 
proton spectral index Tp = 2.0). This energy needs 
to be increased to 1.0 x 10 ^^ (nn/l erg, if 


the cut-off energy in the proton spectrum is reduced to 
30 TeV. When assuming a different proton spectral in¬ 
dex, the required energy in interacting protons increases 
significantly to e.g. Wpp = 25 x 10^^ (nn/lcm”^)”^ 
for the Ep = 2.4 case. 

A source of uncertainty in the leptonic scenario is 
the different radiation fields that contribute to the In¬ 
verse Compton process. One component is the 40 K 
radiation held from 30 Dor C, which we assume to 
have a constant energy density of 0.5eVcm“^. The 
second component originates from the 30 Doradus re¬ 
gion. Depending on the integration region of the Taran¬ 
tula nebula, and the definition of the background re¬ 
gion, the radiation held energy density varies between 
~0.5 eV cm“^ and ~1.5 eV cm“^. Note that these es¬ 
timates assume that the Tarantula nebula and 30 Dor C 
are at the same projected distance. The derived mag¬ 
netic fields for these radiation fields vary between 10/iG 
and 18/iG. The energy in electrons is 3.8 x 10'^® erg to 
2 X 10"^® erg for the different models. The X-ray data 
for 30 Dor C are from (26) and Fetmi-LAT limits on the 
GeV 7 -ray flux have been derived as discussed above. 

The modelling of the N 132D SED has been per¬ 
formed in a similar way to the modelling of 30 Dor C. 
Radio data are from (33), whereas for the upper limit 
on the X-ray synchrotron flux from N 132D we reana¬ 
lyzed Chandra data, and based our estimates only on the 
(4.5 — 6.3) keV band, which is poor in X-ray line emis¬ 
sion (see below). The SNR shell is interacting with a 
nearby molecular cloud, which is also bright in infrared 
wavelengths. Given the complex morphology and lim¬ 
ited angular resolution of HESS, it is not clear where 
the HESS emission is coming from. We therefore as¬ 
sume that all the infrared emission (which has a total 
flux of ~3 Jy, (49)) peaks at about 20 fim (^145 K tem¬ 
perature) and is located in the shell, which has a radius 
of 11 pc. The actual infrared radiation energy density 
experienced by the VHE electrons depends on the loca¬ 
tions of the VHE electrons, which could be the forward, 
but perhaps also the reverse shock, with respect to the 
location from which infrared emission is emitted. If the 
VHE electrons are located all around the forward shock, 
and the IR emission emerges from within N 132D all IR 
photons can in principle be up scattered. Here we as¬ 
sume that about 50% of the photons encounter regions 
with VHE electrons. Under this assumption the lep¬ 
tonic scenario for the combination of 7 -ray flux and the 
X-ray synchrotron flux upper limit requires a magnetic 
field of 2QuG. Ignoring the IR radiation fields, and thus 


http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
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only considering the Cosmic Microwave Background 
Radiation (CMBR), the derived magnetic field is 15/xG. 
The 20/tG field is roughly a factor two lower than the 
equipartition magnetic field in this object. Note, how¬ 
ever, that a detailed modelling of the radiation field is 
beyond the scope of this paper. All input parameters 
for the non-thermal emission model are summarized in 
Tabled 

A.4 X-ray Analysis of 30 Dor C and 
N132D 

The X-ray image of 30 Dor C shown in Fig. 1 in the 
main article was obtained from XMM-Newton obser¬ 
vations (Observation IDs 113020201 and 104660301). 
Only the MOSl and MOS2 detectors were used. The 
data reduction was carried out using the XMM-Newton 
software of the XMM-Newton Science Analysis Soft¬ 
ware (SAS, http://xmm.esac.esa.int/sas/). Calibrated 
event lists were produced for each exposure using the 
SAS emchain script. The periods affected by soft proton 
flaring were excluded (62). To generate images in the 
energy band 0.5 — 8 keV, we created the quiescent par¬ 
ticle background (QPB) images, the count images and 
model exposure maps for each observation and each in¬ 
strument, using mos-spectra and mos back scale. The 
combined mosaic map was finally computed taking into 
account the different efficiencies of the two instruments. 
The image has been smoothed with a Gaussian with a 
width of 10 ". 

The upper limit on the X-ray synchrotron emission 
of N 132D in the 4.5 — 6.3keV continuum dominated 
band is based on the Chandra X-ray observation of Jan¬ 
uary 9, 2006 (ObsID 5532) processed with the CIAO 
v4.5 software. We conservatively assumed that all the 
continuum is X-ray synchrotron emission, whereas in 
reality a thermally dominated origin for the continuum 
emission is more likely given that the spectrum is rich 
in line emission. 

A.5 Radiation Fields in N157B 

The magnetic field estimates presented in the N 157B 
section in the main article sensitively depend on the 
measured 7 -ray and X-ray spectra, as well as the ra¬ 
tio between the energy densities of magnetic fields and 
radiation fields. Depending on the location of N 157B 
in the LMC along the line of sight, different radiation 
fields can potentially act as target for the IC scattering 


Radiation Field 

Temperature 

K 

Energy Density 
eV cm“^ 

CMBR 

2.7 

0.26 

30 Doradus 

88 

2.7 

LH 99 (cold dust) 

29 

12.7 

LH 99 (warm dust) 

230 

5.7 


Table 3: Radiation field energy densities and tempera¬ 
tures. 

of high-energy electrons accelerated inN 157B. Besides 
the omnipresent CMBR, also the infrared radiation from 
the 30 Doradus star-forming region and the OB associa¬ 
tion LH 99 can potentially contribute (if not dominate). 
Table [^summarizes the properties of all radiation fields 
that potentially act as targets for the high-energy elec¬ 
trons accelerated in N 157B. 

We study two extreme cases with maximum and 
minimum radiation field energy densities to derive up¬ 
per and lower limits on the magnetic field in the PWN. 
In Scenario A it is assumed that N 157B is related to 
LH 99 and is located at the same projected distance as 
this OB association, implying that the CMBR, as well 
as radiation associated with 30 Doradus and LH 99, 
contribute to the target radiation fields. In Scenario B, 
N 157B is assumed to be unrelated to LH 99 and to have 
large distances both to the OB association and to 30 Do¬ 
radus so that only the CMBR contributes. 

Several authors argue that N 157B is located in or 
behind the OB association (27, 43, 44), which would 
support Scenario A. In this case the dust emission of 
LH 99 would be the dominant radiation field. The peak 
of the corresponding infrared emission is centred on 
the HII region associated with the young stellar object 
2MASS J05375027-6911071, which is at a projected 
distance of ^ 1' from N 157B. As shown in (63), a two- 
component (cold dust and warm dust) modified black 
body model best describes the spectral energy distri¬ 
bution of the infrared emission from the OB associa¬ 
tion. The LH 99 components have been modeled by the 
modified blackbody model proposed in (63), whereas 
the CMBR and 30 Doradus components are modeled as 
pure blackbody emitters. Since N 157B is potentially 
inside the 2MASS J05375027-6911071 HII region, the 
derived energy densities were furthermore scaled up by 
a factor of 1.5, as would be expected for an object at the 
edge of a homogeneously emitting sphere (64). 

Figure [^ shows the spectral energy distribution of 
N 157B from radio wavelength to 7 -ray energies and 
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Parameter 

age 

yrs 

Eg.p 

E, 

TeV 

B 

pG 

nn 

cm“^ 

VEe 

10 ^® erg 

Wp 

10 ^® erg 

Ti,T2 

K 

Ui, U2 

eV cm“^ 

30 Dor C 

6000 

2.0 

100 

15 

40 

0.38 

15 

88,40 

1.5, 0.5 

N132D 

2350 

2.0 

30 

20 

40 

1.0 

25 

145 

1.0 


Table 2; Input parameters used for the time-dependent modelling of 30 Dor C and N 132D. The electron/proton 
distribution is modelled as dN/dE cx E/Ec). T and u are the temperatures and energy densities 

characterizing the infrared radiation fields. In 30 Dor C two radiation held components are considered, while for 
N 132D, only one component is considered. In all cases, the energy density of the CMBR (1 eV/cm“^) is also 
taken into account for the Inverse Compton process. 



Figure 6 : Spectral energy distribution of N 157B. Also 
shown are the two scenarios as discussed in the text 
for maximum (minimum) radiation helds as straight 
(dashed) lines. A maximum (minimum) magnetic held 
of 45 pG (8 pG) is required to still ht the X-ray and 7 - 
ray data, respectively. For the high-radiation held model 
a power-law index of injected electrons Fg 1 = 2.0 and 
energy in electrons Cep = 0.07 of the current pulsar 
spin-down power is assumed. For the low-radiation 
held model, re ,2 = 2.35 and Ce 2 = 0.45 are as¬ 
sumed. In both scenarios, the cut-off energy is hxed 
to Ee,c = 100 TeV. 


the broadband synchrotron and IC emission expected 
for the time-independent injection of relativistic elec¬ 
trons (400 GeV < Eg < 100 TeV) over the lifetime of 
PSR J0537—6910. In Scenario A an injection of elec¬ 
trons at a constant rate of 7% of the current spin-down 
power of PSR J0537—6910 is required to explain the 
SED, i.e. PF/Enow = 0.07. In this case, a magnetic 
held of 45 pG is required to ht the data. Since the cool¬ 
ing time of electrons that produce the observed X-ray 
and TeV emission would be very short (i.e. < 300 years 
for Eg > 20 TeV), the spin-down power of the pulsar 
should not have changed signihcantly on the scale of 
the electron cooling time. Consequently, this magnetic 
held estimate is rather insensitive to the power-law in¬ 
dex or cut-off energy of the injected electrons, as well 
as to the age or braking index of the pulsar. In contrast, 
in Scenario B, a much smaller magnetic held strength is 
required to explain the data. At the same time, the power 
injected into the nebula needs to be increased to ^45%. 
The smallest held strength still compatible with the X- 
ray and TeV 7 -ray data is 8 p,G, although the spectral 
shape at keV energies is poorly reproduced in the low- 
held case. This mismatch could be overcome if a distri¬ 
bution of magnetic held values is considered, instead of 
a single-magnetic held. The power injected into the neb¬ 
ula of N 157B in this scenario is comparable to the Crab, 
however, the magnetic pressure would be much lower (a 
factor ~240). In this scenario, electron cooling times 
are much longer (i.e. < 8400 years for Eg > 20 TeV) 
and the model is less reliable. 

Note that the modelling of the 7 -ray spectrum is in¬ 
sensitive to the exact choice of the low-energy cut-off in 
the electron spectrum (400 GeV in this case) as H.E.S.S. 
spectral points only start at ^800 GeV. Even higher- 
energy electrons are responsible for the observed X-ray 
emission. The required fraction of pulsar spin down in¬ 
jected into the nebula, however, does change with the 
minimum electron energy. 


15 










A.6 OnthefluxofN157B 

The differential flux of N 157B of $(1 TeV) = (1.3 ± 
0.1) X 10“^^ cm“^ TeV“^) increased by 64% com¬ 
pared to the previously published value (75). A reanal¬ 
ysis of the data set which was used in the previous pub¬ 
lication yields spectral results compatible to the values 
reported here. Thus, the change of the reported flux is 
not intrinsic to the source, the 7 -ray emission of N 157B 
should be considered as being constant. 

The results presented here and in the previous publi¬ 
cation (75) use an image-fitting analysis (18) which pro¬ 
vides an improved angular resolution and a higher 7 -ray 
efficiency. But this reconstruction is also much more 
susceptible to imperfections in the detailed modelling of 
the instrument response than the classical Hillas-based 
analysis. One such imperfection was a misaligned cam¬ 
era which is now corrected in the analysis. In the initial 
publication (75) this correction was not yet taken into 
account. The misaligned camera results not only in an 
underestimated extend of the point-spread function but 
also in a global shift of the ShowerGoodness, the main 
event selection parameter (18). The global shift of the 
ShowerGoodness parameter leads to a misclassification 
of actual 7 -rays as background and thus an underesti¬ 
mation of the true 7 -ray flux. 

This systematic error is relevant only for a small 
number of publications, which utilise the image-fitting 
analysis (18). A re-analysis of affected sources is cur¬ 
rently underway and errata will be published in due 
time. 
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